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We study a periodic one-dimensional exclusion process composed of a driven and a diffusive part. 
In a mesoscopic limit where both dynamics compete we identify bulk-driven phase transitions. We 
employ mean-field theory complemented by Monte-Carlo simulations to characterize the emerging 
non-equilibrium steady states. Monte-Carlo simulations reveal interesting correlation effects that 
we explain phenomenologically. 
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One-dimensional (ID) driven diffusive systems have 
not only served as fruitful testing grounds for fundamen- 
tal questions in non-equilibrium physics Q but have also 
been the focus of recent interest in applications relevant 
to biologically problems Restricted ID motion may 

either result from geometric confinement as in a nuclear 
pore complex of cells |^ or artificial crystalline zeolitical 
structures 0, or arise because molecular engines move 
along one-dimensional tracks as for example in intracel- 
lular transport 1] and protein synthesis . 

In these systems interesting collective effects emerge 
since mutual passage of particles is excluded. The nature 
of these effects depends on vifhether the system is purely 
diffusive {passive) or driven due to the presence of an ex- 
ternal field or an internal driving mechanism inherent to 
the particles (active), e.g. motor activity in intracellular 
transport. In addition, one has to distinguish between 
open and closed boundary conditions. While indepen- 
dent of the boundary conditions both active and passive 
systems show interesting dynamic anomalies [HJ, |ll|, , 
only active systems with open boundaries are known to 
exhibit non-trivial nonequilibrium steady states [l^ . 

Inspired by traffic of molecular motors in closed com- 
partments |3| and colloidal motion in optical traps 
|l2l |. we present a model that combines the symmet- 
ric (SEP) and totally asymmetric simple exclusion pro- 
cess (TASEP). It is intended to investigate the competi- 
tion between driven and diffusive motion in ID systems. 
While the breaking of translational invariance is known 
to be necessary for non-trivial steady states to evolve, we 
show that additionally time scale separation between the 
two processes requires a mesoscopic scaling to guarantee 
a finite current and a physical behavior in the continuum 
limit. 

SEP and TASEP serve as the two paradigms for pas- 
sive and active transport in one dimension. In these 
lattice gas models particles occupy the sites of a one- 
dimensional lattice subject to the simple exclusion rule 
that each site may be occupied by at most one particle. 
In the SEP particles jump independently and randomly 
at rate D to vacant neighboring sites with equal proba- 
bilities to the left and right, while hopping at a rate R 



is strictly unidirectional for the TASEP. In the following 
we will measure time in units of 1/R, i.e. we set R = I. 
For a closed ring geometry both processes are character- 
ized by a steady state with a uniform density profile since 
translational invariance gives equal weight to all permis- 
sible configurations. For open boundary conditions, the 
steady state of the SEP always shows a linear density 
profile , whose slope depends on the boundary condi- 
tion's difference. In stark contrast, TASEP exhibits sev- 
eral distinct nonequilibrium phases [l^ as a function of 
the magnitude of the entrance and exit rate, a and /3, at 
the left and right boundary, respectively. Two different 
phases can be characterized by their global particle den- 
sity (low and high) and one by a maximal current. The 
validity of the initial mean-field analysis of the phase dia- 
gram was later supported by several exact solutions |l5l | . 

In this Letter we aim at identifying the nature of 
the nonequilibrium steady states of a closed ring sys- 
tem consistiirg of two equally sized compartments with 
i = I,-- - ,N lattice sites, whose dynamics is governed 
by a symmetric and a totally asymmetric exclusion pro- 
cess, respectively; see Fig^ To distinguish between the 




FIG. 1: Schematic model of the ring system. The dynamics 
of the upper (passive) lane with lattice sites Xi is governed by 
a SEP with rate D, and the dynamics of the lower (active) 
lane with lattice sites Xi by a TASEP with unity jump rate. 

active and passive compartments of the ring, quantities 
for the passive part like the location of the lattice sites, 
Xi, or the occupation numbers hi S {0,1} are indicated 
by a tilde. Particle exchange between both sub-lattices is 
exclusively allowed at their junctions with dynamic rules 
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defined by the originating site. Thus a vacancy on the 
right hand side ii of the passive part can be filled with 
a particle from site xjv of the active section with rate 
unity, while the corresponding event at the left junction 
will occur with rate D. Since the system is closed, the 
particle number A^p = + "-0 conserved and 

the particle density 
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serves as a dimensionless control parameter confined to 
the interval [0, 1] . Having chosen the inverse hopping rate 
1/R on the active part as our time unit, we are left with 
a two-dimensional parameter space (rip, D). As these pa- 
rameters are bulk quantities, resulting phase transitions 
will be bulk-induced in contrast to boundary-induced 
phase transitions in TASEP. In order to explore the 
phase behavior in this parameter space we will exploit 
the fact that for a closed system the steady state current 
J{np,D) is spatially constant. 

We analyze the ring system by mean-field (MF) the- 
ory complemented by Monte Carlo simulations (MCS). 
The MF analysis allows to decouple the two parts of the 
ring and consider them as separate lanes with effective 
entrance and exit rates. Once these rates are identified 
we can use known results for the SEP and TASEP with 
open boundaries. To begin with, we introduce the fol- 
lowing notation for the stationary densities pi = (rii) at 
the junction sites: 7 = pi, S = ppf for the passive part, 
and a = pi and 1 — /? — pn according to TASEP con- 
vention for the active part. The incoming current to site 
xi is the product of the average occupation number on 
the originating site, the jump rate to the destination site 
and the probability that the latter is empty: 5D{1 — a). 
Due to current conservation this has to equal the current 
to site X2'. a{l — a). Using particle- hole symmetry, one 
can proceed accordingly at the other junction to arrive 
at the following relations: 



a^DS, 7 = 1-/3 
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To analyze the interplay of the two junctions, we will 
exploit the conservation of current mentioned above. On 
the active part, directed motion results in a current of 
J = pi(l — Pi+i). Conservation of current can only 
be fulfilled for a spatially constant or piecewise con- 
stant density distribution. In the latter case two sec- 
tions of constant density U and V are connected by a 
domain wall and their densities have to fulfill the condi- 
tion pu = \ — pv to conserve the current. 

As particle motion is bi-directional on the passive part, 
the current between two sites is obtained as the balance 
of their bilateral particle exchange which is proportional 
to their density difference. Conservation of current thus 
demands a linear density slope and the current takes the 
form J = {j — S)D/N reminiscent of Pick's law Evi- 
dently the passive part current vanishes with system size 
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FIG. 2; Simulated density profiles of the active part in a sys- 
tem of A*' = 500 and Tip = 0.38 for different diffusion strengths 
d (indicated in the graph) . The system is in the LD-HD phase 
featuring a domain wall. With increasing d the domain wall 
position (intersection with the solid line) is shifted to the left 
(see Eq. ||KJ) and DW fluctuations increase. The left bound- 
ary layers are due to current limitation caused by correlations 
at the passive part boundaries (see text). 



while the active part current is constant. Due to cur- 
rent conservation this implies that the smaller current 
and thus the diffusive process is dominating the system 
in a trivial way. To broaden our analysis to a wide pa- 
rameter range for arbitrary system sizes we introduce a 
mesoscopic scaling ^ and a new control parameter: 
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This scaling can be understood as a time scale separation 
and ensures competitive behavior between the system's 
constituents. The specific form Q even guarantees a 
well-behaved approach to the continuous case N —f 00. 
With the new control parameter the passive part current 
is expressed as J = d{-^ — d) and allows for a relation 
between the two junctions: 
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Having derived suitable entry and exit rates for the ac- 
tive part, we can now apply the TASEP results. The low 
(high) density phases (LD, HD) are realized in the peri- 
odic system equally and are characterized by a uniform 
density below (above) 1/2 and a boundary layer at the 
right (left) . Special attention has to be paid to the phase 
boundary a = (3 between LD and HD phase, where the 
boundaries are matched by a piecewise constant density 
profile with an intervening domain wall (DW). Due to 
the randomness of entry and exit events this DW is dclo- 
calized and subjected to a random walk that explores the 
complete system on long time scales. In the ring system 
however, our MCS reveal a DW localization (Fig|21l. We 
understand this as a consequence of entry and exit rates 
that are not statistically independent as in TASEP but 
rather connected via the passive part, dependent on the 
amount of available particles and the diffusion rate. 
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To proceed with a quantitative analysis, we connect 
the boundary conditions of the two sub-lattices in a self- 
consistent manner to obtain a junction density solely de- 
pendent on the two control parameters. The precondi- 
tion a = /3 guarantees a conserved active part current 
of J = a{l — a). This allows to rewrite Eq. I^J as 
6 = (1 — a)[l — a/d]. Since S is of order 1/N in the 
TL, we find a = d + d 0{1/N). Having derived a, the 
remaining densities f3, 7 and S are easily calculated and 
we can proceed to connect the position of the DW to 
the global particle density. To this end continuous den- 
sity distributions are assumed which are applicable for 
large systems where the lattice spacing vanishes. To be- 
gin with, a DW is assumed to be present on the active 
part. Hence, a Heaviside function that connects two re- 
gions of constant density at DW position is chosen as 
p{x) — a + Q{x — Xvj){l — a — (3). For the passive part's 
density distribution a linear slope p{x) = (5 -f (7 — 6)x is 
appropriate. Now particle conservation can be expressed 
as 2np = dx[p{x) -\- p{x)]. Solving this in the TL the 
DW position evaluates to 

_ -3 + 3d + 4np 

The DW position is left to depend on the two control pa- 
rameters: the diffusion strength (see Fig. |21for an exam- 
ple) and the particle density. The resulting (d, np)-phase 
diagram exhibits several phase transitions and regimes 
of which three can be characterized by the domain wall 
position. In the case x^ <Q the DW has left the active 
part at the left junction resulting in an high density phase 
(HD) with constant density p{x) = 1-/3. A constant low 
density (LD) corresponding to the left boundary condi- 
tion p{x) = a is estabhshed for 1 < Xw, while the DW 
is localized inside the system for < Xw < 1 connect- 
ing the two boundaries by a phase of coexistence (LD- 
HD). Its phase boundaries can be obtained from Eq. Q 
as d = (3 — 4np)/3 for Xw = and d — 4rip — 1 for 
Xw = 1 • The phase boundaries intersect at a critical 
point at (rip — 3/8, d = 1/2) where the DW height van- 
ishes (Fig. 13). The extremal points at d = of the 
LD-HD phase can be readily explained by the particle 
fraction of a quarter that has to be bound on the pas- 
sive part for vanishing diffusion (signifies a = /3 = 0). 
Hence, at rip = 1/4 there are no particles available to 
the active part and the system crosses over into the LD 
phase. Analogous, the HD phase is entered for values of 
rip > 3/4. The existence of the LD-HD phase is a dis- 
tinctive difference to TASEP . It replaces one discon- 
tinuous phase transition with two continuous transitions 
that meet at a multi-critical point. 

Similar to TASEP, the ring system also features a 
maximal-current phase (MC). In this regime the active 
part current imposes its maximum Jmc — 1/4 on both 
lanes. MF approximations allow again to derive the ex- 
tent of the MC phase in dependence of the two control 
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FIG. 3: Phase diagram obtained by MF (solid lines) exhibits 
four phases. Simulations (MC (triangles) and LD-HD (cir- 
cles) phase boundaries, N = 200) reveal a failure of the MF 
analysis: correlations at the passive part boundaries cause a 
diminished current that shifts the MC phase. 

parameters d and rip. To this end we have to assume a 
different density distribution than above for the active 
part. The corresponding distribution in TASEP is es- 
tablished for boundary conditions a,/3 > 1/2 and has 
to be constant at p = 1/2 with the exception of possi- 
ble boundary layers that vanish in the TL. Hence, A^/2 
particles have to be present on the active part. We can 
then deduce from particle conservation: 2np — 1/2 — 
(5 + (7 — (5)/2 where the r.h.s is just the integral over 
the passive part's linear density distribution. Using the 
equality of passive part and active part current allows 
to solve for 5 — 2 (rip — 1/4) — l/(8d). The active part 
constraints can be rewritten as constraints on the passive 
part by use of Eq. Q and Eq. (0) to 5 > l/(2diV) and 
5 < 1/2 — l/(4d). The last three equations constitute 
planes in the ((5, d, np)-space. Computing now the inter- 
section of the first plane with the two inequalities, one 
can finally deduce the phase boundaries of the MC phase 
in (d, np)-phase space in the TL as d = l/(16np — 4) and 
d — 1/(8 — 16np) where the latter is the boundary with 
the HD phase. Notice that the former phase boundary 
could already be derived by considering only terms of 
order 1 in the TL (i.e. with 5 = 0) while the HD bound- 
ary is obtained only US — 0(1/ N) is considered. The 
MC phase originates at the critical point in phase space 
and asymptotes for d ^ 00 at values of 1/2 and 1/4 (see 
Fig. as can easily be explained. At particle densities 
below Hp — 1/4 there are not enough particles available 
to establish a constant density of p = 1/2 on the active 
part. On the contrary, the system makes a transition to 
the HD phase if both active and passive part are half- 
filled at Up — 1/2. At this point 7 > 1/2 implies a 
violation of the MC phase requirement (3 > 1/2. 

We have complemented our MF analysis by extensive 
MCS. Their results deviate from the phase diagram de- 
rived above (see simulation data in Fig.OJ. The key dif- 
ference is that the multi-critical point is shifted to higher 
values of d. To rationalize this we compare the density 
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profiles from simulations and MF theory. While the den- 
sity distribution on the passive part has been assumed 
to be linear according to SEP results, simulated data ex- 
hibit a distinct curvature in the profile [2^. A closer 
examination of the present boundary conditions reveals 
that these are different from open boundaries used for the 
derivation of the linear SEP density profile. As particles 
enter the passive part with a rate of unity while the in- 
ternal dynamics of the diffusion are much faster, the SEP 
experiences a boundary that is refiective during a consid- 
erable amount of time. Like in other problems, exhibiting 
time scale separation, as forest-fire models [l^ the sta- 
tionary state results from competition between a process 
that strives towards equilibrium (diffusion) and another 
process repeatedly driving the system to nonequilibrium 
(particle entrance). The rare particle entrance events 
provide a particle excess at the right boundary. This 
density surplus then spreads into the bulk by diffusion 
before the next entrance restores the particle excess [l7j . 
A time average then results in a strictly convex density 
profile on the passive part. Its curvature renders Pick's 
Law for the passive part current inappropriate because 
the smallest gradient (here at the left interface) acts as a 
current bottleneck. The consequences for the phase dia- 
gram are obvious: since the current is effectively reduced 
the MC phase is only established at higher diffusion val- 
ues (Fig. El- The current hmitation caused by the bot- 
tleneck is also responsible for the boundary layers on the 
l.h.s. of the active part (Fig|21). Here global conservation 
of current forces the density to decay quickly to a value 
conform with the system's minimal current. Thereby the 
extension of the LD-HD phase to values of c? > 1/2 can 
be explained. By introducing a correction factor for the 
reduced current a refined MF theory can explain the shift 
of the MC phase at least qualitatively. Recalculation of 
the LD-HD phase boundaries in this realm completely 
fails, also due to correlation effects on the left interface. 

Our analytical and numerical studies have shown that 
the presented system exhibits a rich phase behavior un- 
expected for periodic systems. While boundary-induced 
phase transitions in exclusion processes occur in several 
geometries 0, ll^ and non-trivial nonequilibrium steady 
states are known in periodic systems if translational sym- 
metry is broken by defects ll8| , comparable bulk- induced 
phase transitions have not been studied to our knowl- 
edge so far. Crucial to this behavior is the competition 
between the two processes which the scaling Eq. © en- 
sures for a broad parameter range. The time scale sep- 
aration between the subprocesses is also responsible for 
the quantitative failure of MF that comes as quite a sur- 
prise, since up to now similar approximations have been 
known as a reliable tool to reproduce the phase diagram 
of lattice gas systems with an astonishing accuracy. 

For the system presented, these characteristics indi- 
cate the existence of interesting correlation phenomena 
that call for analytical methods beyond MF. The phe- 



nomenological explanation and the modified MF expres- 
sion for the current should therefore only be considered 
preliminary. As MCS near the multi-critical point re- 
quire considerable computation resources, we cannot give 
any reliable predictions about the exact phase topology. 
The exact form of the phase diagram and in particu- 
lar the behavior near the multi-critical point, remains 
an open question to whose solution different techniques 
like Bethe ansatz or density matrix RG may 
contribute. Furthermore, possible experimental realiza- 
tions [1^ could make this system an intriguing problem 
to study. 
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